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Qh! Abstract. 2-point topological charge correlation functions of several types of 

' geometric singularity in gaussian random fields are calculated explicitly, using a general 

t:;^ . scheme: zeros of n-dimensional random vectors, signed by the sign of their jacobian 

I determinant; critical points (gradient zeros) of real scalars in two dimensions signed 

' by the hessian; and umbilic points of real scalars in two dimensions, signed by their 

^ , index. The functions in each case depend on the underlying spatial correlation function 

' of the field. These topological charge correlation functions are found to obey the first 

■ Stillinger-Lovett sum rule for ionic fiuids. 
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^ ' 1. Introduction 

Although a spatially extended field may be smooth, and contain no infinities, it may 
. nonetheless have point singularities associated with its topology. For example, a 2- 
^ dimensional landscape, specified by a real function, has critical points (stationary 

points), where the gradient of the field vanishes and the gradient direction cannot be 
defined; in the neighbourhood of such a point, the gradient direction changes arbitrarily 
quickly through all of its values. Such points are characterised by a topological charge, 
a signed integer which is determined by the local geometry of the singularity; in the 
case of critical points, the number is the (signed) number of rotations of the gradient 
vector around the critical point; it is +1 for maxima and minima, and —1 for saddles. 
The topological charge is also called the Poincare index (or Hopf index) (Milnor 1965), 
and is defined as a signed integer at the point zeros of n-dimensional real vector fields 
in n dimensions. Although the number of zeros may change as a field evolves, the total 
topological charge is constant; it is a topological invariant. 

Here, I discuss the statistical properties of various zeros/singularities, in fields which 
are specified by gaussian random functions. In this case, only structurally stable zeros 
of charge ±1 occur. Specifically, I shall discuss how densities and two point charge 
correlation functions of the distributions of signed singularities may be calculated under 
a rather general scheme, and then use this scheme to calculate the topological charge 
correlation functions of three types of singularity in isotropic random fields: zeros of n- 
dimensional vectors in ra- dimensional space (rederiving a result originally due to Halperin 
(1981)); critical points of random scalar fields in two dimensions; and umbilic points 
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of random scalar fields in two dimensions. The charge correlation functions for each 
of these types of topological singularity are found to be different in each case, and 
dependent on the underlying correlation function of the gaussian field. Each function is 
found to satisfy a screening relation associated with ionic liquids. 

Topological zeros are very important in many areas of physics and mathematics: 
in addition to critical points (gradient zeros) which have obvious importance, the zeros 
of 2-dimensional complex scalar fields (phase singularities, wave dislocations, vortices, 
realised as 2-dimensional vectors) are also of great interest, especially where the field 
is a quantum wavefunction or an optical field (e.g. Nye and Berry 1974, Berry and 
Dennis 2000, Dennis 2001b), or an order parameter (Mermin 1979). Only point zeros 
are considered here - the dimensionality n of the vector field, whose zeros provide the 
singularity, is assumed equal to the dimensionality of the space. 

The first systematic study of the statistical geometry of random real scalar fields 
in two dimensions was by Longuet-Higgins (1957a,b, 1958), who generalised one- 
dimensional methods of Rice (1954) to calculate, amongst other things, the densities 
of critical points and probability density function of the gaussian curvature of the 
function. Halperin (1981) derived the n-dimensional vector correlation function, whose 
proof was supplied by Liu and Mazenko (1992), and has recently been recast in the 
language of riemannian geometry by Foltin (2003a). Various correlation functions in two 
dimensions, including those for phase singularities and critical points, were investigated 
numerically by Freund and Wilkinson (1998). Planar phase singularity correlations 
(including density correlations) were investigated by Berry and Dennis (2000). 

The topological singularities of interest are the zeros of an n-dimensional vector field 
V = {vi, . . . ,fn). The field is a smooth function defined on an n-dimensional euclidean 
space, with points labelled by the vector r = (ri, . . . ,r„). In random fields, the only 
statistically significant zeros are those of first order, whose jacobian determinant 



is nonzero (where di = •^i = 8/ dr.)). The topological charge of such zeros is given by 
When the field is random, the density d{r) of zeros at a position r is 



where (•) denotes averaging over the statistical ensemble, and the zeros are identified by 
the n-dimensional Dirac 5-function. The modulus of the jacobian is included to ensure 
that each zero has the correct statistical weight, so d{r) has the units of density. This 
expression confirms that degenerate zeros are so rare they have no statistical significance. 
The average topological charge q{r) at r is expressed as the density ()1.2|) . but each 
singularity is weighted by its charge, that is, the sign of the jacobian. The charge 
correlation functions g with which this paper is concerned are the generalisation of the 
average charge density to that at two points r^i, r^, normalised by the density: 



J = det diVj, 






9{rA,rB) 
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In the following section, the scheme for calculating topological charge correlation 
functions ()1.3|) in general gaussian random fields is explained. Section El then provides 
background to gaussian fields which are statistically stationary and isotropic. The 
scheme is then applied to zeros of n-dimensional vector fields (section H}, critical points 
of 2-dimensional scalar fields (section IH)), and umbilic points of the same fields (section 
in} ■ The phenomenon of topological charge screening in these three cases is discussed in 
section [7| 

As this paper was being completed, I became aware of the work of Foltin (2003b), 
who performs similar calculations for critical points by a different, possibly simpler 
method. 

2. Gaussian random functions and scheme for calculating charge 
correlation functions 

An ensemble of scalar functions / in n-dimensional space is said to be a gaussian random 
function if the probability distribution of the function at each point of space is given by 
a gaussian distribution, so the function at each point defines a gaussian random variable 
(Adler 1981). The only restrictions on gaussian random functions made in this section 
are that they be centred (the average (/) = 0) and that first derivatives exist. 

Our starting point is the well-known expression for the probability density function 
for a set of independent gaussian random functions u = {ui, . . . ,un), 



It is now possible to present the general scheme for calculating 2-point topological 
charge correlation functions of zeros of gaussian random functions, expressed in (|1.3|) . 
The topological charges considered are all zeros of a real vector gaussian random function 
V = (fi,...,f„), of dimension n, and the correlation function is the average of the 
product of the local density at points labelled A, B. Quantities evaluated at these places 
are denoted with the appropriate subscript A oi B. The variables that appear in the 
average ()1.3j) are therefore the components of va.vb, and their derivatives Ova^Ovb 
that appear in the jacobians J{va) = J a-, J{,vb) = Jb- There are a total of m different 
derivatives diVj appearing in each jacobian, and in general m < n^. For example, in 
the case of critical points, n = 2, and m = 3 (the terms appearing in the jacobian 
are c^n/, ^22/, ^12/ = 821 f). The calculation therefore involves an average of the 
A^ = 2(m + n) -dimensional gaussian random vector 



The correlation matrix S for this u is defined as in ()2.2|) . and averages are evaluated 
according to the probability density function (j2IH). 




(2.1) 



where S is the correlation matrix, with components defined 



(2.2) 



U = {{dVA)l, {dVA)m, {dVB)l, {dVB)m, Vai, • • • , VAn, Vbi, • • • , Vbu) ■ 
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The correlation function in (jl.3|) is therefore 
{6^{vaW{vb)JaJb) 
d{rA)d{rB) 



d{rA)d{rB) (27r)"+'"VdetS 

X j d2("+™)ii5"K)5"(i^B)JAjrBexp(-ii ■ ■ u/2). (2.4) 

where (i(r^), d{rB) are the appropriate densities of zeros ()1.2p . The jacobians Ja, Jb can 
be quite comphcated; each is a muhihnear function, involving a sum of n-fold products 
of Ml, ... , Um (for A) or Um+i, • • • , U2m (for B). 

The 2n x 2r;,-dimensional lower right submatrix of S (i.e. the averages dependent 
only on VaiVb) is denoted by K; the complementary 2m x 2m-dimensional upper left 
submatrix of is denoted by H~^; that is, 

(2.5) 



K 

where the other blocks marked • are not denoted by special symbols. In [Appendix A[ 
Jacobi's determinant theorem is used to show that 

det S = det K det H. (2.6) 

The integrals in vai-, ■ ■ ■ , vbh, only involving 5-functions, are performed, leaving the 
2m derivative terms dv to integrate; let u' = {{dvA)i, • • • , {dvB)m)- Then 

9AB = I <^''''u JaJb^M-u' ■ ■ u' /2) 

ci(rA)rf(T'B)(27r)"+'"VdetS J 

1 



d{rA)d{rB) (27r)"+™Vdet S 
V det H 



X 



j d^-^u' JaJb^^^^ j d^™* exp(it ■ W - t ■ H • t/2), (2.7) 

upon Fourier transforming the gaussian in u' with the 2m-dimensional Fourier vector 
variable t. The jacobian terms Ja (depending on u'^, . . . ,u'^) and J7b (depending on 
M^^^, . . . , M2m,) may be replaced with partial derivative terms in tj, 

u'j iVt^. = iV,-, (2.8) 
where V is used to denote partial derivatives in 2m-dimensional £-space. The jacobians, 
with this replacement, have now become differential operators, denoted Jva,J\7b- 
Therefore 

1 /detH 
^^"^ ~ rf(rA)c/(rB)(27r)-+2-V d^ 



X 



J d^'^t exp(-t ■ S • t/2)J^AJvB J d'™tx' exp(i£ ■ u') 

1 (27r)2™ 



d(r^)c/(rB)(27r)"+2- 



X 



J d^^t exp(-i ■ H • t/2)J^AJvB^^'^{t), (2.9) 
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where, in the second hne, the integral over u' is reahsed as the Fourier transform of 
the (5-function of t, and ()2.6|) has been used in the prefactor. The expression is then 
integrated by parts (so the Jsj operators act on the exponential term rather than the 
^-function), and then integrated in t. The final expression is 

= d{r^)dirs){2.rV^^' ^'-'^^ 

where 

D = {J^aJ^b exp(-i ■ H • V2)]t=o • (2-11) 
The charge correlation function therefore only depends on d, det K and components of 
the inverse reduced inverse correlation matrix 3. 

This final step, evaluating D in ()2.11|) is the most complicated part of the 
calculation, and depends on the precise form of the jacobian determinant J . Each 
summand in the operator JsjaJvb is a 2n-fold derivative over some index set {r} = 
{ti, . . . , T2n}i and is comprised of n terms from Vi, . . . , Vm (for and n terms from 
Vm+i,---,V2m (for B), either set possibly including repetitions. The result of this 
particular operation is 



2n 



WiS/r, exp(-t ■ H ■ t/2) 



all pairings 



= E n ^-^v' (2.12) 
t=o P (-fA,TA/)ep 

where the sum on the right hand side is over all pairings p = (tcTq,') . . . (r^,r^') of 
indices in {r}; there are {2n — 1)!! such pairings, each involving n sets of pairs. The 
n-fold product on the right hand side is over all components of H with indices given by 
the appropriate pair. The Hjj themselves are found by further application of Jacobi's 
determinant theorem in [Appendix A , and are expressed in terms of minors of det S in 
(ESI). 

The correlation functions calculated explicitly in this paper are not too complicated, 
either because n is small (only two-dimensional fields are considered in sections El); 
or H is sparse, as in the case of random vectors (section |3]) . 

The density d of zeros appears in ()2.10p : in general, this may be difficult to calculate, 
due to the modulus sign in ()1.2j) . For random n-dimensional vector fields, the main 
part of the density calculation is in [Appendix B[ For critical and umbilic points, these 
densities were calculated by Longuet-Higgins (1957a,b) and Berry and Hannay (1977). 



3. Isotropic random fields 

The topological charge correlation function in ()2.10j) is extremely general, applying to 
any centred differentiable gaussian random vector field. In this section, and for the 
remainder of the article, attention will be restricted to stationary isotropic random 
fields. For these fields, all averages are (statistically) translation and rotation invariant. 
They are conveniently given by a Fourier representation 

/(r) = 5^a(A;)cos(fc-r + 0fc), (3.1) 

k 
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where k are now the Fourier variable vectors (wavevectors) . The components Vi of 
random vectors are specified by independent identically distributed realisations of ()3.1|) . 
The amplitude a{k) only depends on the magnitude |A;| = k,s and the phase (pk is 
uniformly random - each ensemble member is therefore labelled by the choice of (pk for 
each k. It also may represent the spatial part of a real linear homogeneous nondispersive 
wavefield, for which the representation ()3.1|) is particularly evocative. The infinite k set 
is assumed sufficiently dense that they may be represented as an integral, and 

J2(^^ik)»^ f (TkUik); (3.2) 
fc 

where Il{k) is the power spectrum of the field; by the Wiener- Khinchine theorem (Feller 
1950), n(fc) is the n-dimensional Fourier transform of the field correlation function C(r), 
where r = — r^l, and 

C{r) = {fAfs) = {f{rA)f{rB)), (3.3) 

normalised such that C(0) = (/^) = 1. The only condition on C is that it is symmetric 
and has positive Fourier transform. 

Averages of derivatives of / may be represented as moments of 11 (as by Longuet- 
Higgins 1957a,b, Berry and Hannay 1977, Berry and Dennis 2000), or equivalently 
derivatives of C as is done here. Coordinates are chosen where 

r = rB- Va, ri = r, = 0, j = 2, . . . , n. (3.4) 

Since the fields are isotropic, the results are not affected by this choice. The correlations 
computed in sections El are in two dimensions; in this case, direction 1 is denoted by 
X, direction 2hy y. 

Representing derivatives by subscripts, the correlations of first derivatives of / are 
found to be (Berry and Dennis 2000) 

E ^ {fAAfs) = -(/a/b,i) = - = -C 

H ^ ifAjB,,) = -{fA,nfB) = -{fAfB,n) = - [dlC]^^^^^^^^^^^^^ = -C'/r, j ^ 1. (3.5) 

Averages involving / and its first derivatives other than those in ()3.3p . ()3.5p are zero. The 
averages are equal to appropriate derivatives of C(r), and then setting ri = r, r2, . . . = 
(as in ()3.4j] ). and derivatives in of odd order vanish. The functions in ()3.5j) . when 
r = (denoted by subscript 0), are 

= 0, Fo = H^ = -C"(0) = -Cl (3.6) 

It is easily verified that —Cq > 0, since C has a positive Fourier transform. The 
correlation matrices ()5.4|) . ()6.5p required to calculate the topological charge correlation 
functions of critical points and umbilic points, involve higher derivatives of C, given in 

(lEsi), dnsi). 

As the separation between A and B increases, the correlation between fA and 
decreases and C ^ 0. This decay is slowest in the case where all the k in ()3.1|) have the 
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same length k^, (set to 1 for convenience) and the power spectrum Il.{k) = 6{k — 1). For 
any n > 2, the corresponding correlation function is a Bessel function times a factor of 
r, with oscillatory decay that falls off like r~*^"~^)/^. Of particular interest is the n = 2 
case, for which C(r) = Jo{r). The spectrum in this case was called the ring spectrum by 
Longuet-Higgins (1957a,b) and Berry and Dennis (2000), and is conjectured to model 
the high eigenfunctions in quantum chaotic billiards (Berry 1978, 2002). 



4. Correlations of zeros of isotropic vector fields 

In the present section, we shall consider n-dimensional gaussian random vector fields 
V = (fi,...,f„) in n dimensions whose cartesian components are independent and 
identically distributed gaussian random fields (the derivatives of the components are 
also assumed completely independent). The jacobian J', whose sign determines the 
topological charges of the zeros, is the determinant of the matrix of first derivatives 

We begin by calculating the density p.2|l of zeros of random vectors. This was 
calculated for = 1,2,3 by Halperin (1981) and Liu and Mazenko (1992), (the n = 1 
case was previously found by Rice (1954), and n = 2 hj Berry (1978)). For general n, 
the density dn ()1.2|) is expressed as a probability integral with density function 1)2.11) 
and correlations given by ()3.5|) . ()3.6p . Therefore 

dn = {6''{vi)\detVij\) 



1 



(2vr)"("+i)/2Fo" 

_ ^0 



/ 1 " 1 



(27r)"("+i)/2 




2 



(T'v,,, I det Vij I exp - - V vl. , (4. 1) 



where in the final line the (5-functions in the Vi were integrated, and the Vij were rescaled 
(each by y/Fo) to be dimensionless. The remaining integral is solved in [Appendix B[ 
giving the density of zeros in n dimensions 

^ _ f n"\nl2 ( 2 )• _ / r^n\n/2 ~ l)!o"n-l „n 

In this expression, cr„_i is the surface area of the unit {n — l)-sphere in n dimensions, 
given by ()B.4|) . As is common in such problems in statistical geometry, the result is a 
spectral quantity ((— C'q)"/^) multiplied by a geometric factor. 

The scheme of the previous section may now be applied to calculate the 
topological charge correlation function for zeros in these gaussian random vector fields. 
Implementation of the scheme is facilitated by the fact that the components Vi are 
completely independent, and m = n^. The submatrix K„ of the correlation matrix S„ 
only depends on the correlations of the components of the vectors va,vb', from section 
[HI the only correlations that do not vanish are {v\j) = (f^j) = 1 and {vaiVbi) = C for 
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2 = 1, . . . , n. It is easy to see that 

detK„ = (l-C^)". (4.3) 

From ()2.12|) . the other necessary ingredient of the correlation function scheme is the 
components of the matrix S, defined in ()2.5j) . The elements of 3 are labelled by the 
multiindex of the components VAi,j,VBk,h using the correlations ()3.5|) . ()3.6|) and the 
arguments of Appendix A[ particularly equation ()A.3|) . the only nonvanishing elements 
are 



^{Ai,l){Ai,l) — ^{Bi,l)(Bi,l) 
^{Ai,l){Bi,l) = ^{Bi,l){Ai,l) 



'^{Ai,j)iAi,j) - '^{Bi,j){Bi,j) - Fq 

'^(Ai,j){Bi,j) = '^{Bi,j){Ai,j) = H for i = 1, . . . , n, j = 2, . . . , n. (4.4) 
The problem remains to use these components and ()2.12j) to evaluate ()2.1H) . By 

(EH), (El), 

JvA = det VAij, (4.5) 
(similarly for JTvb)- Each of the n\ summands in this determinant is an n-fold product 
signer ^ Ai,a{i)) where a is a permutation of 1, . . . , n. Therefore 

permutations n 

JvaJ^b= ^ signer signcr'(-l)" JJ VAi,a(i)VBj,a(i)- (4-6) 
(t,(t' *,i=i 

From ()2.12|) . the result of one of these summands acting on exp(— t ■ St/2) and setting 

t = is nonzero if there is a pairing of these multiindices where the corresponding 

elements of H are nonzero. From ()4.4j) . this is only the case when the permutations a, a' 

are the same. Thus, from ()2.12|) and ()4.4j) . 

permutations n 

Yl Yl-iAi,amBi,aii)) 
a i=\ 



'T-! W '^{Ai,i){Bi,i) 



i=l 



= n\{F -CE^/{l-C^))H'^-^. (4.7) 
This, together with ()4.3|) . can now be put into ()2.10p to give 

n! (F(l - C^) - CE^) [ H ^ 



9n[r) 



2ndl (l-C2)3/2 \2'kVT^ 
n\ (C"(l - C^) + CC'^) [ -C" 



27rd2 (1 - C2)3/2 \2TxrVT^ 
in-l)\ dhJr) 



n-l 



(27r)"d2^"-i dr 
where the function /?.„(r) is defined 



(4.8) 
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a b 




Figure 1. The 2-diniensional vector zero correlation function g2{r), plotted for two 
choices of C(r) : (a) C(r) = Jo(r); (b) C{r) ^ exp(-r2/2). 

This decays to as r — oo, and when r = 0, 

KiO) = i-C'^r/' (4.10) 

(the quantity Cq = C"{0) is always negative, since C is the Fourier transform of the 
positive power spectrum 11). 

The topological charge correlation function of zeros of n-dimensional gaussian 
random vector fields g{r) was first calculated by Halperin (1981), in a form equivalent 
to ()4.8p . but without proof. This function was also derived by different means by Liu 
and Mazenko (1992), and in the n = 1 case by Rice (1954), and the n = 2 case by 
Berry and Dennis (2000) and Foltin (2003a). g{r) is plotted in figure Qfor two choices 
of the field correlation function C for n = 2. When C{r) = Jo{r), g2 is oscillatory; when 
C(r) = exp(— r^/2), it increases monotonically to zero as r — oo. 

5. Critical points in two dimensions 

In this section, the scheme of section E] is used to calculate the topological charge 
correlation function of critical points of isotropic gaussian random functions in the plane, 
that is the Poincare index correlation function of isotropic random surfaces. 

The gaussian random function examined shall be written / = f{r) = f{x, y) (where 
X = ri,y = r2)i which will be assumed stationary and isotropic, so the expressions in 
section El may be used. In particular, / and its first derivatives have the correlations 
fl3.5|) . ()3.6|) . as well as further correlations involving second derivatives, described below. 
As in the previous section, for convenience in calculations, the two points A and B are 
separated only in the x coordinate. 

At a critical point, the gradient V/ = (/,x, f,y) is zero. The critical point jacobian 
Jc whose sign defines the topological charge (Poincare index) is the hessian determinant 

J7c = det f^ij = f^xxf,yy ~ f,xy'j i^-^) 

this is the gaussian curvature of the surface. Unlike a general 2-dimensional random 
vector field (as in the previous section), the gradient field V/ is irrotational, which gives 
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rise to relationships and correlations between the derivatives of the components (e.g. 
f,xy = f,yx^ whereas before, fi_2 7^ f2,i in general). This makes the explicit computation 
of the topological charge correlation function more difficult than in the previous section, 
even in the 2-dimensional case that is considered here; the scheme of section |21 applies 
to gradient zeros in fields of any dimension. 

The statistical properties of critical points of gaussian random fields in 2 dimensions 
were considered by Longuet-Higgins (1957a,b); he found that the density dc of critical 
points is (Longuet-Higgins 1957b equations (71), (78)): 



dr 



2C, 



(4) 



(5.2) 



(Cq'''* denotes the fourth derivative of C, evaluated at r = 0). The density of saddles 
equals the density of extrema (maxima and minima), and the density of maxima equals 
the density of minima. The probability density function of the gaussian curvature J , 
despite its asymmetry (Longuet-Higgins 1958 equation (7.14), Dennis 2002 equation 
(57)), has zero first moment, implying the average topological charge f)Jc) is 

zero, as expected. 

The topological charge correlation function is again calculated using the scheme of 
section|2l particularly (j2.10p . Therefore, the vector u (12. 3p of gaussian random variables, 
in a convenient ordering, is 



nc = {fA,xx, fA,yy, fB,xx, fB,yy, fA,xy, fB,xy, fA,x, fB,x, fA,y, fB,y) 



(5.3) 



rix (c. 
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(5.4) 



where the correlations between elements of are computed to be 



I ^[dldyG]^^^^^^, 
Go = lo = 0, 



C(3) 

(rG" 



2G')/r' 



=r,y=0 
=r,y=0 



= 3{rG" 
Mo = No 



G')/r^ 
3Lo = 



G, 



(4) 

• 



(5.5) 
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The last line gives the special value of these correlations when r = 0. 

The matrix is the 4x4 lower right submatrix of Sc? and has determinant 

detK, = (F2-i72)(^2_^2)_ ^5Q) 

The pair of differential jacobian operators are, from ()5.1|) . 

JvaJvb = V1V2V3V4 + V^V^ - ViVsV^ - V3V4V^. (5.7) 

Using ()2.12|1 . the result of these operators acting is on exp(— t ■ S ■ t/2) and setting 
t = (c.f. (imUl . dTTTD ) is 

Dc = '='14'='23+'='13'='24~2^16^26 + '='12'='34 — 2^35^45— ^34^55 + 2^50 — ^12^66+'='55'='66; (5-8) 

where the necessary entries of 3c, are found using Jacobi's determinant theorem in 
l73|) : as an example. 



-24 



Mo-FP/{F^ -F^). (5.9) 



The topological charge correlation function gdr) for critical points is obtained 
by substituting (Ol) (with all terms like found using (I A. 3D into (ICTHl ). 

This expression is complicated and not very illuminating, and is not given here. Upon 
substituting ()3.5|) . ()3.6p . ()5.5p in, one finds that gc can be written as a perfect derivative 

(c.f. 63), 

where 

iC" - C'/r) 

hc{r) - 



C(3)(3Cf -2C"2-C"C7r 



jC" - CVr)(3Cf - 2C^Vr^ - C"C'/r) 

r{C[;^ - cyr^) 

1 d iC'-C'/rf 

r{C" - C'/r) dr _ C"2)(c^'2 _ ^,2/^2) " 

This process of finding he is long and tedious, and details are omitted here. It is easy 
to see that /ic ^ as r ^ cxd; it is straightforward, by Taylor expanding derivatives of 
C, to show that 

4(7(4) 

hc{0) = = -2nd.. (5.12) 

The critical point topological charge correlation function for two particular field 
correlation functions is shown in figure |21 As with the 2-dimensional vector case plotted 
in figure Q the properties of the correlation function (oscillatory, exponential decay, 
etc) are similar to that of the underlying field correlation function C(r), on which the 
correlation function depends. Features of interest in these plots include the sharp initial 
minimum of (a), and the fact that the monotonic decay in (b) is from above, not below 
(in contrast to its counterpart in figure ^ Nevertheless, the form of h. is significantly 
more complicated than especially when n = 2. 
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a 





Figure 2. The critical point charge correlation function gdr), plotted for two choices 
of C(r) : (a) C{r) = Jo(r); (b) C(r) = exp(-rV2). 



6. Umbilic points 



Less well-known than critical points, umbilic points are geometric point singularity 
features associated with the second derivative of / - namely where hessian matrix of 
second derivatives dijf is degenerate (Berry and Hannay 1977, Porteous 2001, Hilbert 
and Cohn-Vossen 1952). Geometrically, the principal axes of gaussian curvature are not 
defined at these points. The eigenvalues of the hessian coincide when f^^x = f,yy, f,xy = 
0. Umbilic points are therefore zeros of the 2-dimensional vector field 

Vn= iif,xx- f,yy)/2J,^y). (6.1) 

The factor of half in the first term ensures that is statistically rotation invariant. 

An umbilic point has an index, determined geometrically by the sense of rotation 
of the principal axes of curvature around the umbilic point, and the index is generically 
±1/2 (Berry and Hannay 1977); only the sign of the index is of interest here, and this 
is determined by the appropriate jacobian J7u on Vu, 

2<^u = f,xxxf,xyy + f,yyyf,xxy ~ f,xyy ~ f ,xxyi (6-2) 

depending on the third partial derivatives of /. The calculation of the topological charge 
correlation function for umbilic points can proceed according to the scheme of section 
13 in a similar way to the corresponding calculation for critical points. 

Umbilic points for isotropic random functions was considered by Berry and Hannay; 
they found the density d^ to be (Berry and Hannay 1977 equation (34)) 



3C, 



(6) 



IOttC, 



(4) 



(6.3) 



(Cq^'* is the sixth derivative of C at 0) and that the average index q is zero (the separate 
densities are O.Sdu for stars, 0.05279(iu for monstars, and 0.44721(iu for lemons). In the 
present work, the distinction between monstars and lemons, which both have positive 
index, is not used. 
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The ordering of the vector of gaussian random functions liu is chosen 

= fA,.yy. fs,..., fB,.yy, fA,..y, fA,yyy, fB,..y, fB,yyy, 

{fA,xx ~ fA,yy)/'^y {fB,xx ~ fB,yy)/'^y f A,xy, fB,xy)- 



(6.4) 



The correlation matrix ()2.2j) is 
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Po = Qo = Ro = 0, 




5*0 — 5To — 


5f/o 


= K) 













(6.5) 



The matrix Ku, defined using p.5|l . has determinant 

det K^ = {Ll-{M + N- 2T) V16)(T2 - L^) 
The result of the jacobian derivative operators 1)2.111) gives 



(6.6) 
(6.7) 



T, = H 



12 



^2 
"14 



^2 
"22 



2"24 ~ 2^12^22 + '='13"24 



— 2.^oo^S 



4^14^24 — 2^12^55 + 2^12^56 
-22'='56 + "55 + "56 + 2^57 + ^53 — 2^55^56 — 4^57^53 + ^57^68- (6.8) 

The necessary entries of Hu are found using ()A.3|) . The resulting expression for D^, and 
therefore for gu_, is very complicated, but may be reduced to the following form: 

1 dhu 



where 



9n[r) 



hJr) 



Air'^d'^r dr 



d_ f r{Q-Rf 
dr V 4 Vdet 



Q{P + R-2Q) 
4 Vdet Ku 



(6.9) 



(6.10) 
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a b 




Figure 3. The umbilic charge correlation function gu{r)^ plotted for two choices of 
C{r) : (a) C(r) = Mr); (b) C{r) = exp(-rV2). 

It can be shown that /iu(0) = — 27r(iu, and h ^ Q asr ^ oo. gn{f) is plotted in figureElfor 
C(r) = Joir) and exp(— r^/2). The most obvious feature of these two plots, compared to 
figures n 121 is that they have a negative maximum, a property that seems to be general 
for (7u; although this has not been proved. It is unclear what the physical significance 
of this kink should be; mathematically, it probably arises from interference between the 
two summands in h^{r) in ()6.10|) . 



7. Topological charge screening 

Three particular charge correlation functions have been derived exactly (equations ()4.8|) . 
(I4.9|l . (jSmj, In each case, the charge correlation function has the 

form 

gir) = ^ (7 1) 

where n is the number of dimensions, d is the density of zeros and h{r) is a function 
such that 

{2'KYd 

/i(0) = — ^—^ -, /i oo asr ^ oo. (7.2) 

(T„_i(n - 1)! 

The total charge density around a given (say positive) topological charge is therefore 
d / dV^(r) = V^^n-i / 



(27r)"rf ""Vo dr 
in - l)!o-„-i 



-(/i(oo) - /i(0)) 



(27r)"d2 

= - 1, (7.3) 

where the hypersphere area ct„_i appears in the first line from integration in polar 
coordinates. It implies that the distribution of topological charges is such that every 
topological charge tends to be surrounded by charges of the opposite sign, such that the 
topological charge is 'screened' at large distances. This fact was noticed in the random 



Correlations and screening of topological charges in gaussian random fields 



15 



vector case by Halperin (1981) and Liu and Mazenko (1992) (although the general zero 
density dn had not been determined explicitly) and is independent of the field correlation 
function C. The derivation here shows that this is a more general phenomenon, possibly 
a universal feature of topological charge correlations for gaussian random fields. It 
should be noted that ()7.Hj] is not necessarily satisfied for an arbitrary distribution of 
signed points; for instance, g{r) = always for Poisson points, for which there is no 
screening. 5-function correlations at the origin are ignored in the following. 

An analogy may be drawn from the theory of ionic liquids (Hansen and McDonald 
1986); in a fiuid or plasma, consisting of two species identical except for opposite 
(Coulomb) charges, the following Stillinger-Lovett sum rules (Stillinger and Lovett 
1968a,b) are found to hold: 



djd-rgir) =-1, (7.4) 

d j dVr2^(r) = -a„Al (7.5) 

Here, g is the charge-charge correlation function, A is a characteristic screening length 
dependent on temperature and ra, and constant dependent on dimensionality. 

These rules are discussed for n = 3 by Hansen and McDonald (1986), Stillinger and 
Lovett 1968a,b, and n = 2 by Jancovici 1987, Jancovici et al 1994. 

The screening relation ()7.3j) is equivalent to the first Stillinger-Lovett sum rule ()7.4|1 . 
which is derived using the statistical mechanics of pairwise. Coulomb interacting fiuids. 
It is unclear whether the fact that topological and coulombic charges screen in the same 
way is coincidence, or evidence of some deeper connection between the two statistical 
theories. 

It is natural to ask whether the topological charge correlation functions satisfy the 
second sum rule, which (upon integrating the left hand side of ()7.5p by parts), depends 
on the integral of rh. For n = 2, it was found (Berry and Dennis 2000) that for certain 
choices of C, this integral may diverge. The slowest decay comes when C = Jq, and by 

h2{j') ~r-^oo cos(r + 7r/4)^/r when C = Jq, (7.6) 

giving a logarithmic divergence for the second moment. For comparison, the critical and 
umbilic functions hch^, (equations (jSHH), ()6.1()j) respectively) both give, for the same 
choice of C, 

/ic(r), /iu(r) ~r^oo cos(2r)/r^ when C = Jq, (7.7) 

implying that the second sum rule is satisfied generally for critical and umbilic points, 
although, as for the cases of random vectors where the integral converges, the screening 
length A, defined in analogy to ()7.5j) . depends on the choice of C. For random vectors 
with n > 2, the second moment of g always converges, because of the higher power of C 
appearing in ()4.9|1 (also, the decay of C will be faster, as discussed at the end of section 
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Comparison may be drawn to the electrostatic analogy in random matrix theory, 
particularly in the case of the so-called Ginibre ensemble of NxN matrices whose entries 
are independent, identically distributed circular gaussian random variables (Ginibre 
1965). The eigenvalues of these matrices are found to have exactly the same statistical 
behaviour as a 1-component 2-dimensional Coulomb gas of N charges in a harmonic 
oscillator potential, and the 2-point density correlation function screens against a 
uniform background (i.e. d J d^2; ((^cin ^ 1) = ^1); and have finite second moment. 
Certain random polynomial analogues have zeros that can also be expressed as 2- 
dimensional Coulomb gases with additional interactions (Hannay 1998, Forrester and 
Honner 1999). The eigenvalues of random matrices (which may be expressed as the 
zeros of the characteristic polynomial) and zeros of random polynomials are all of the 
same sign, since they are zeros of complex analytic functions, and the density correlation 
functions of zeros are unique (there is no analogue of C). 

There may be a danger in taking the analogy with fluids too far; for instance, the 
oscillations of the functions in figures lip2pHb are reminiscent of those of charged fluids 
(e.g. Hansen and McDonald 1986). However, the physical causes for these oscillations 
are very different. In fluids, the oscillations usually arise from packing considerations 
(the ions themselves are of finite size, fixing the lengthscale, although in plasmas they 
are usually represented as point charges (Bans and Hansen 1980)). Topological charges, 
on the other hand, are points, and the oscillations in these figures originate from the 
oscillations in the underlying field correlation function C(r), which is Jq in this case. 
Although Halperin (1981) discusses the similarity between the short-range behaviour of 
the 2-dimensional vector correlation function ()4.8|) . ()4.9|) and Kosterlitz-Thouless theory, 
the present situation is more general, both in that the (possibly long-range) screening 
is exact, and that the results hold for any reasonable field correlation function C{r). 

8. Discussion 

Using a general scheme for calculating topological charge correlation functions, three 
particular correlation functions were found explicitly, and were found always to satisfy 
a screening relation ()7.3|) . 

The scheme of section |21 used to calculate the charge correlation functions is very 
general, and can be generalised to calculate the charge-charge correlation function 
between two different sets of topological charges - for instance with a critical point 
at A, and an umbilic point at B. Although not done so here, the scheme may be applied 
to anisotropic fields. 

It is rather more difficult to calculate the density correlation function for topological 
charges (the analogue of (jl.Hj) where the moduli of the jacobians are taken). It was 
calculated for 1-dimensional vectors (i.e. random functions in 1 dimension) by Rice 
(1954), and for 2-dimensional vectors (realised as complex scalars) by Berry and Dennis 
(2000), Saichev et al (2001). However, it has not been possible to generalise such 
methods to the case of critical points. Also, numerical evidence (Freund and Wilkinson 
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1998) suggests that the correlation function of extrema signed by the sign of their 
laplacian (+1 for minima, —1 for maxima) also satisfies the screening relation ()7.4j) . All 
such functions would be needed to calculate the partial correlation functions between 
the species (e.g. maxima with maxima, maxima with saddles, maxima with minima, 
etc), which would give a complete statistical picture. 

The presence of boundaries in the random function will affect the statistical 
properties of topological charges (e.g. Berry (2002) for nodal points in the plane), 
and it is possible that there may be some further analogy with the physics of interfaces 
of Coulomb fluids. The scheme of section |21 ought to be adaptable to calculate charge 
correlation functions in this case. 

Only zeros of fields linear in gaussian random functions have been considered here, 
although the density of others may be calculated, for instance, in addition to nodal 
points, a 2-dimensional gaussian random complex scalar has critical points of its modulus 
squared (Weinrib and Halperin 1982) and its argument (Dennis 2001a). The scheme 
employed here cannot be used to calculate correlation functions for these, although 
numerical evidence (I Freund, personal communication) suggests that the critical points 
of argument (together with the nodal points) do screen, therefore adding to the cases 
shown here. It is tempting to conjecture that topological charge screening may be a 
universal phenomenon in gaussian random fields. 
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Appendix A. Jacobi's determinant theorem 

Let A be a square matrix. The minor 7Viy^'"*^(A) is the determinant of the k x k 

submatrix of A with rows ii, . . . ,ik and columns ji, . . . ,jk- A^*^^' ( A) shall be used 
to denote the complementary minor, that is, the minor of the submatrix of A with rows 
ii, . . . ,ik and columns ji, . . . ,jn excluded. Then Jacobi's determinant theorem (Jeffreys 
and Jeffreys 1956, page 135) states 



Applying this to S in ()2.5|1 . and choosing H ^ as the submatrix whose determinant is 
the minor of 



detAA<-::-(A- 




(A.l) 



det S det H 



1 



detSA^l;:::;^™(S-^) 




, ,^2m+l,...,2{m+n)/-sp\ 
■'^^2m+l,...,2(m+n)\^) 



detK, 



(A.2) 



from which ()2.6p follows directly. 
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Jacobi's theorem can also be used to find the elements Sjj of the inverse reduced 
inverse matrix H in ()2.5|) . needed for ()2.12|) . In this case, ()A.1|) is applied twice, once 
on the matrix H~^, and once on S. Therefore 

= (-l)^+^A7j.(S-^)/detH-i 

= (-l)^+^detH^i;:;P\^;;;:|:,(S-^) 

det" A^i.2m+l,...,2m+2n/y,N 
uei ^ J^'- j,2m+l,...,2m+2n\^ ) 



« yi,2m+l,...,2m+2?i) 
•'^j,2m+l,.--,2m+2n 




= detK-Met y 'J ^ j (A.3) 
where in the last line • represents the terms Sj fc, j where k = 2m + 1, . . . , 2m + 2n. 



Thus the Sjj appearing in the expression for D in ()2.12|1 . is the determinant of the 
(2n + 1) X [2n + 1) submatrix comprised of the ith row and jth column of S, and the 
submatrix K. 



Appendix B. Calculation of the density of vector zeros (jHH) 

In order to integrate 1^4.1^ . the following must be integrated 

V = I d"%,, I dett;,,| exp (^-^ ■ (B.l) 

This is, mathematically, the average hypervolume of an n-dimensional parallelepiped 
specified by gaussian random vectors wi = {vi^i, . . . vi^n), ■ ■ ■ ,Wn = . . . , f^.n)- 

These gaussian random vectors are identically distributed isotropically in n-dimensional 
space. The hypervolume is nonzero exactly when the set of n vectors is linearly 
independent. 

This hypervolume may be found explicitly in a manner reminiscent of the Gram- 
Schmidt orthogonalization procedure for vectors. Geometrically, 

volume of parallelepiped = length of Wi x length of W2 orthogonal towi x . . . 

X length of lUn orthogonal to span{ to 1, . . . ,iUn-i}- (B-2) 

A given factor in this product is therefore the average length of the gaussian 
random vector Wj in the orthogonal complement of a (j — l)-dimensional subspace 
spanjtui, . . . , Wj^i}. 

Now, since the vector Wj is isotropic, it may be represented identically in any choice 
of orthonormal basis of n-dimensional space; in particular, its first j — 1 components 
t>j 1, . . . , t>j j_i may be chosen to be in sp3Ji{wi, . . . ,Wj_i} (as in the Gram-Schmidt 
procedure). The total contribution of Wj to the integral in ()4.H1 involves the average 
length of the vector made up of the other components Vjj, . . . , Vj^n- Where k = n—j — 1, 
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this is 



oo 1 

-oo 



d'-'^k-i / dpp'=exp(-p72)(B.3) 



where, in the second hne, the first m — 1 components have been pulled out as trivial 
gaussians, integrating to (2txY^~^^/'^ = (27r)'^"~'^)/^; the remaining k integrals are the 
average length of a gaussian random vector in fc-dimensional space. This integral has 
been converted to polar coordinates, with dtUfc-i the solid angle infinitesimal on the unit 
{k — l)-sphere Sk-i, and p is the radius. It is well known that the surface area crfc_i of 
the unit [k — l)-hypersphere is 

d'-'n.^, = (B.4) 



The p integral in ()B.3|) is 2'^'^^^^/^ ((A; — l)/2)!. Therefore, the numerical part of ()4.H) is 
l/(27r)"'*^""'"^)/^ times the product in ()B.2|) . with each term in the product, now labelled 
by k, given by the expression ()B.3|) . Therefore 



V 



(^2Ti)n(n+l)/2 



(27r)"("+l)/2 

(¥)! 



" 9„A:/2 

J](27r)("-'=)/2 X 4^ X 



k=l 



)/2^n/2 



k=l 



(¥)' 
(¥)' 



2(A:-l)/2 



k-\ 



(B.5) 



^(n+l)/2 ■ 

This value agrees with that stated by Halperin (1981), Liu and Mazenko (1992) of l/vr 
(n = 1), l/27r {n = 2) and l/vr^ (n = 3). 
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